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Abstract. In this paper we report the theoretical results obtained for partial ionization yields and the 
above-threshold ionization (ATI) spectra of Magnesium in a Tksapphire laser field (804 nm) in the range 
of short pulse duration (20-120 fs). Ionization yield, with linearly polarized light for a 120 fs laser pulse, 
is obtained as a function of the peak intensity motivated by recent experimental data [9]- For this, we 
have solved the time-dependent Schrodinger equation nonperturbatively on a basis of discretized states 
obtained with two different methods; one with the two-electron wavefunction relaxed at the boundaries, 
giving a quadratic discretized basis and the other with the two-electron wavefunction expanded in terms 
of Mg + -orbitals plus one free electron allowing the handling of multiple continua (open channels). Results, 
obtained with the two methods, are compared and advantages and disadvantages of the open-channel 
method are discussed. 

PACS. 3 2.80.Rm Multiphoton ionization and excitation to highly excited states (e.g., Rydberg states) 



1 Introduction 

Experimental studies of above threshold ionization have 
concentrated mostly on the rare gases [l][2l[3j|4] which can 
withstand sufficiently large intensities of radiation of wave- 
length around 780 to 800 nm, which is the most versatile 
and convenient source of short pulse duration. As a re- 
sult, extensive studies of ATI photoelectron spectra, and 
of course high order harmonics, in the non-perturbative 
regime have and are being conducted [5,6,7,8. Softer atoms, 
with lower ionization potentials are not expected to ex- 
hibit the extensive spectra of ATI, as they ionize at rela- 
tively lower intensities. Nevertheless, recent experimental 
studies [MTU] have addressed this issue in atomic Mag- 
nesium, which not only is a soft atom but also has two 
valence electrons, which may or may not add additional 
features, depending on intensity and pulse duration. The 
results have been rather surprising. The slopes of the ATI 
spectra as a function of laser intensity are difficult to un- 
derstand p. This motivated the theoretical study whose 
main results are presented in this paper. The study was 
undertaken because we happened to have developed the 
theoretical tools necessary for a realistic description of 
Magnesium as a correlated two-electron system. Although 
under the conditions of the experiment [5] , it does not ap- 
pear that two-electron excitations played an important 
role, at least in the aspects we have addressed, still we 
can say that the atomic structure entering the calculation 
is realistic. We have not calculated everything that has 



been measured experimentally. We present our results as 
a point of calibration of what can be expected theoreti- 
cally. As discussed in the sections that follow, we find it 
difficult to understand the observed slopes. If it is because 
something is missing in our calculation, at this point we 
are unable to say what it could be. 



2 Two-electron continuum states of 
Magnesium 

The Hamiltonian of the Magnesium atom H a can be writ- 
ten as: 



-^ + K//(r i ) 



+ V(r 1 ,r 2 ), (1) 
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where V e ff (rj represents the effective potential for the i th 
electron and the core, and V(ri, T2) is a two-body interac- 
tion operator, that includes the static Coulomb interaction 
1 / 1 r i — r2 1 as well as a dielectronic effective interaction [11] 

UJ. 

The characteristic feature of Magnesium is the exis- 
tence of an ns 2 valence shell, outside a closed-shell core, 
the excitation of which requires a much larger amount 
of energy compared with the first and second ionization 
threshold. This allows us to explore excitation and/or ion- 
ization processes of the valence electrons, without consid- 
ering the closed-shell core-excitation, for a certain range 
of photon energies. The construction of the Magnesium 
bound and final states is achieved through a standard 
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Hartree-Fock calculation and the inclusion of a core-po- 
larization potential which represents the influence of the 
core on the two valence electrons in a way very similar to 
that described in (TTJEI]. 

At a first stage, we perform a Hartree-Fock calcula- 
tion for the closed-shell core of Magnesium (Mg ++ ), deriv- 
ing thus the effective Hartree-Fock potentials 'seen' by an 
outer electron. At a second stage, we solve the Schrodinger 
equation for Mg + . In this case, the effective potential act- 
ing on the valence electron is given by [12] . 



V ef f(r) = Vr(r) + 



1 



-(-) E 



(2) 



where a s is the static polarizability of the doubly-ionized 
Mg and r; cut-off radii for the various partial waves I = 
0,1,2,.... 

Having produced the Mg + one-electron radial eigen- 
states P n i(r) for each partial wave I = 0,1, 2, we solve 
the Schrodinger equation: 



H a # A (r 1 ,r 2 )=EV A (r 1 ,r 2 ) 



(3) 



by expanding the two-electron eigenstates ^ /1 (ri,r 2 ) on 
the basis of two-electron antisymmetrized orbitals namely: 

^ A (r 1 ,r 2 ;E i ) = £ Cmvn'iEi)^,^,^) (4) 

nil 1 ,n' 



';'( r i, r 2j 



A, 



n r 2 
x Y L M L (r 1 ,r2;l,l')Y s ,M s (s,s'), (5) 



where Ai 2 is the antisymmetrization operator which en- 
sures that the total wave function, is anti-symmetric with 
respect to interchange of the space and spin coordinates 
of the two electrons. A represents the set of angular quan- 
tum numbers (LSM^Ms). In Eq.© the choice of the ba- 
sis functions R n i(r) determines the method used for the 
calculation of the correlated two-electron states. 

In the case of fixed boundary conditions (FXBC), we 
force the wavefunction to be zero at the boundaries by 
selecting the basis functions R n i (r) to be the one-electron 
radial solutions of Mg + , which by construction vanish at 
the boundaries. This, transforms the Schrodinger Eq. (J3J) 
into a generalized eigenvalue matrix equation, by the diag- 
onalization of which we obtain the coefficients C n u> n '(Ei) 
for each discrete eigenvalue Ei [14llT51[T5] . This choice of 
the basis functions R n i (r) leads thus to a discretized con- 
tinuum spectrum for the Magnesium atom with density of 
states determined basically from the box radius R. 

On the other hand, we can select the basis function 
R n i(r) to be non- vanishing at the boundaries (FRBC) 
transforming the Schrodinger Eq.Q into a system of al- 
gebraic equations for the coefficients C n u> n i(Ei) [T51IT5] . 
For this, we use for both cases (FXBC and FRBC) the B- 
spline basis functions 17 used extensively in the last 15 
years in atomic structure calculations. The main advan- 
tage of this method over the FXBC method is that the 



density of the continuum spectrum is completely control- 
lable (therefore degeneracy into the continuum is guar- 
anteed) and the two-electron continuum states have been 
expanded in terms of core-target states plus an ionized 
electron, as in standard close-coupling scattering theory. 

We have ascertained that the results obtained with the 
free boundary method compare well with the convergent 
results obtained with the fixed boundary method. 



3 Time-dependent Schrodinger Equation 

Our objective is to solve the time-dependent Schrodinger 
equation (TDSE) : 

^(n, r 2 ; t) = -i [H a + V(t)] !?(n, r 2 ; t). (6) 

The time-dependent interaction V(t) of the two-valence 
electron atom with an external laser pulse in the dipole 
approximation and velocity gauge can be written as: 

V(t) = -A(t)e-( Pl +p 2 ), (7) 

where the vector potential A(t) is given by, 

A(t) = A f(t)cos(ujt), (8) 

e is the unit polarization vector of the laser field, and pi i2 
are the momenta of the two valence electrons. Aq = Eq/uj 
is the amplitude of the vector potential, and Eo the elec- 
tric field strength. The pulse shape envelope is f(t) = 
cos 2 (irt/r), where r is approximately the full width at 
half maximum (FWHM) for the electric field. The integra- 
tion time for the cosine squared shape pulse, is taken from 
— t/2 to t/2. The velocity form of the interaction Hamilto- 
nian is chosen, because it makes the calculation converge 
faster, in terms of the number of angular momenta in- 
cluded [12] ■ It should be noted here that length and veloc- 
ity dipole matrix elements agree very well. For the length 
operator, for the reason that the core-polarization poten- 
tial is angular momentum dependent, we have checked 
that the use of the modified operator as given by |18j . 
r' = r — f(a s /r 2 )[l — e~( r / r ^ ] has very small effect in the 
dipole matrix elements compared with those calculated 
with the standard length operator. 

The time-dependent wavefunction is now expanded on 
the basis of cigenfunctions x P a {ti, r 2 ; Ei), 



!?(ri,r 2 ;t) = ^ ^(^^(ri.rajEi), 

Ei,L 



(9) 



transforming the time-dependent Schrodinger equation into 
a set of first-order differential equations for the time-dependent 
coefficients Z7 i (£^ i ,f), namely, 

.j t U L (Ei,t) = J] [Ei,LkvSLL'-V Ei L t B i ,L'{t)]UL'{E V ,t) 



En L' 



(10) 

subject to the initial conditions \Ul{Ei, t = 0)| 2 = <5e 4 _Ei<^lo- 
The quantities Ve^.e -,L'(t) represent the matrix elements 
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of V(t) calculated between the states characterized by the 
quantum numbers EiL and Ei'L' . Assuming that the laser 
field is linearly polarized, and the Magnesium atom ini- 
tially in its ground state (L = 0, h = h = 0, M = 0), we 
need to consider only the M = singlet states {S — 0). 

At the end of the pulse, we obtain the coefficients 
UL{Ei,t — > oo), from which we can extract information 
about the ionization yields, photoelectron energy spec- 
trum as well as photoelectron angular distributions [TS1 

EQ|. 



4 Results and discussion 

Compared to the Helium atom, the Magnesium continua 
associated with more than one ionization threshold lie 
much closer to the lowest one (one open channel) which 
makes them much more accessible at intensities, frequen- 
cies, and pulse durations currently available. That is why 
excited ionic states have often been observed in alkaline 
earth atoms. 

For laser intensities of about Io < 10 12 W/cm 2 of a 
Thsapphire laser at 804 nm (1.54 eV), five photons are 
required to ionize the Magnesium atom 

Mg(3s 2 ) + bu -> Mg+(3s) + e"(0.06 eV), (11) 

and leave the Mg + in its ground state, as illustrated in the 
energy-level diagram in Figure [TJ At higher intensities as 
the ponderomotive shift increases, at least six photons are 
required to ionize the Magnesium atom and the 3s3p state 
can be reached through a three-photon absorption. If we 
consider, for example, 11-photon ionization in Magnesium, 
five of the photons will bring the atom in the ground ionic 
state plus one electron in the associated continuum, and 
the remaining ones will be absorbed with the possibility 
to excite further the ionized electron (above-threshold ion- 
ization) and/or to excite the core to different ionic states. 
The calculation of the probability of this channeling into 
the final ionization stages of the Magnesium, after the 
interaction with the laser, is a difficult task. In the follow- 
ing paragraphs, we discuss the results obtained with two 
different methods each having its own advantages and dis- 
advantages. 

Populating the first excited ionic state Mg + (3p) is pos- 
sible by the absorption of eight photons. With the ab- 
sorption of an additional three photons, the ionic excited 
states Mg + (4s) and Mg + (3d) become also energetically 
accessible. The position in energy of the ionic thresholds 
in Figure [U is for the bare atom without taking into ac- 
count the shift of the energy levels in the laser field. In 
the presence of an intense laser, the ionization threshold 
of an atom is increased by the ponderomotive energy. If 
the intensity of the laser is sufficiently high, channel clos- 
ing will occur so that NHuj < I p + U p , where I p = —E a 
is the ionization potential corresponding to the energy of 
the ionization threshold E a , N is the number of photons 
required to ionize the atom, and U p is the ponderomotive 
energy which is given in atomic units by U p = Eg/4w 2 . In 



the case of photon energy ui = 1.54 eV the ponderomotive 
shift is Up ~ 0.6 eV at the intensity 7=10 13 W/cm 2 . 

A related phenomenon which occurs in strong laser 
fields is that large AC-Stark shifts in atomic states can 
lead to intermediate resonances in the ionization process 
[2"T] . At intensities around 2.4 x 10 12 W/cm 2 , the excited 
state of the Magnesium atom 3s3p is shifted into reso- 
nance. 



4.1 Fixed boundary condition method 

We have studied the intensity dependence of both ion yield 
and the photoelectron spectrum of Magnesium in the Ti- 
Sapphire laser field in the domain of short pulse duration 
for laser intensities up to 10 14 W/cm 2 with two-electron 
states of Magnesium calculated with the fixed boundary 
method. The results presented are converged in terms of 
box size, number of B-splines and angular momenta. 

In Figure we show the yield of Magnesium ions as a 
function of the laser intensity, for photon energies u) — 1.54 
eV, and laser pulse duration r = 35 optical cycles (~ 93 
fs) and r = 45 optical cycles (~ 120 fs). We have employed 
a basis of 992 B-splines in a box of 1900 a.u, and a total 
number of angular momenta up to L = 9. 

According to perturbation theory, the A^-photon gen- 
eralized cross section is proportional to the N th power 
of the photon flux. Therefore the slope of the ionization 
yield versus laser intensity in a logarithmic scale repre- 
sents the minimum number of photons necessary to ionize 
the atom. Since the Magnesium states in the continuum 
are shifted in the laser field with the ponderomotive po- 
tential, at least six photons must be absorbed to ionize 
the Magnesium atom for laser intensities higher than 10 12 
W/cm 2 . 

For intensities in the range between 10 12 W/cm 2 and 
10 13 W/cm 2 , and for both the 93 and 120 fs pulses we 
found slopes for the ionization yield curves around 6 which 
are basically in agreement with the lowest-order pertur- 
bation theory (LOPT) predictions. The ionization yield 
of Mg + increases linearly with the peak intensity of the 
field until the saturation regime is reached, where a de- 
crease in the slope begins. Around 10 13 W/cm 2 , however, 
we have found a sudden enhancement of the ionization 
yield which is attributed to AC-stark shifting of the Ry- 
dberg states into resonance with the laser field. Though 
this enhancement is apparent in our results, the slope does 
not compare to that (of about 11) reported in the recent 
experiment by [5], in the same range of intensities. 

In Figure [3] we have plotted the ground state popu- 
lation, the ionization yield as well as the population of 
states 3s3p and 3s4/. In this figure one can see that at 
peak intensity of about 1.3 x 10 13 W/cm 2 the ionization 
yield makes an abrupt 'jump'. At the same intensity the 
population of the state 3s4/ is almost 10% of the ioniza- 
tion yield (almost equal with the remaining population in 
the ground state) while the population of the 3s3p state 
becomes negligible. In Figure [4] we have plotted for the 
peak intensity 1.3 x 10 13 W/cm 2 the populations of the 
same states as in Figure [3] and the ionization yield during 
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as a function of the time. The 3s3p state is populated dur- 
ing the rise of the pulse while at the end of the pulse the 
population of the 3s3p state has been 'transferred' to the 
3s4/ state, which becomes the dominant, in population, 
excited state. 



4.2 Free boundary condition method 

We produce the initial ground state of Magnesium in a box 
of radius R = 40 a.u. in which we have included only 20 
configuration series of the type nsms, npmp, ndmd, nfmf 
with n, m up to 5. The number of B-splines basis was 62 
and the order 9. The parameters of the core-polarization 
potential included in EqJ5]were a s — 0.491 for the static 
polarizability of the closed-shell core and 1.2, 1.335, 1.25, 
1.3, 1.1 for the cut-off radii for the I — 0,1,2,3,4, par- 
tial waves, respectively. This resulted to a value for the 
singly ionized Magnesium equal -E(Mg + (3s)) = —14.99 
eV and a value for the Magnesium ground state equal to 
—22.653 eV measured from the double ionization thresh- 
old. The energy value for the intermediate 3s3p state was 
-18.32 eV. 

The basis used for the propagation of the TDSE differs 
from the one used in the FXBC case in two aspects: firstly, 
we have produced our basis in a much smaller box (of ra- 
dius 40 a.u) than that of the FXBC basis (~ 600 2000 
a.u.) which results in the absence of high Rydberg states 
in the Mg bound spectrum. Thus, while in the FXBC ba- 
sis a number of Rydberg states of about ~ 30 -j- 40 might 
be obtained, depending on the value of the box radius, in 
the FRBC basis a number of 5 Rydberg states is obtained. 
Moreover, we have calculated a total number of angular 
momenta, up to five. The only reason for this restriction 
of the FRBC basis is computational resources limit. It is 
the price that one has to pay in going from the FXBC 
method to the FRBC method. The expansion of the two- 
electron states in a set of target - core states plus one free 
electron, allows us to obtain much more detailed informa- 
tion for the states of the atom after its interaction with 
the laser, such as partial ionization yields (populations of 
the various ionization stages of the ionized system) , par- 
tial photoelectron spectra (branching ratios) and moreover 
angular photoelectron distributions. Thus it is expected, 
that quantitative differences are to be found between the 
two approaches due to the limited basis used in the FRBC 
approach. It is absolutely necessary, therefore, that one 
should ensure that results calculated with both methods 
are in good agreement, since in an unlimited basis case 
they have to be the same. 

With this in mind, we have compared in Figure [5] the 
ionization yield in Magnesium calculated with the free 
boundary method for a pulse duration of about 29.5 fs 
(full curve) with the ionization yield calculated with the 
fixed boundary method of about 27 fs (dotted curve). We 
have employed a basis of 802 B-splines in a box of 900 a.u 
and calculated two-electron states with total angular mo- 
menta up to L = 9, for the results obtained with the fixed 
boundary method. For short pulse duration, the agree- 
ment between these two methods is rather good. The cor- 



responding slope of the ionization yields is 4.46 and 4.25 
for the FXBC and FRBC method, respectively. Similar 
agreement is obtained for the photoelectron energy spec- 
trum for the pulses employed here, though not presented. 

In Figure El we have plotted the ion yield versus inten- 
sity of the laser field, of Magnesium ions left in different 
states, for photon energy w = 1.54 eV and 93 fs pulse 
duration, calculated with a cosine squared pulse. From 
the ionization yield curves, we found that the Mg ion is 
preferentially left in the 3s ground state in the intensity 
domain smaller than 2 x 10 13 W/cm 2 . The slope of the cor- 
responding curve (represented by the dotted line) has the 
value 5.96 around the peak intensity Iq = 10 12 W/cm 2 . We 
should mention here the difference in the slope of the ion- 
ization yield between the shorter pulse (fig. [5]) which gives 
a slope of about ~ 4.4. For the longer pulse the photon en- 
ergy spectrum of the field is much less broader than in the 
case of the shorter pulse. Due to bandwidth of the shorter 
pulse, a three photon absorption from the ground state 
populates strongly the 3s3p excited state which for that 
reason, plays the role of an intermediate resonance thus 
lowering the order of the slope of the ionization yield, as 
expected from energy conservation rules, according LOPT 
predictions. Note that, ionization from the 3s3p state re- 
quires the absorption of three more 1.55 eV photons. Sim- 
ilar role plays the ionization threshold at about ~ 7.64 eV, 
while the 5-photon (1.55 eV) absorption brings the system 
in an energy very close to the ionization threshold. 

In the presence of resonances, the ionization rate for 
multiphoton processes measured as a function of intensity 
will deviate from the power-law I N dependence [22] . Since 
we are dealing with near resonant multiphoton ionization 
with the 3s3p state of Magnesium, the slopes of the ion- 
ization yields for the Mg+(3p), Mg+(4s), and Mg+(3d) 
thresholds are smaller than those predicted by perturba- 
tion theory. In order to have a better picture of the amount 
of ions left in different excited states, we have calculated 
the branching ratios in photoelectron spectra correspond- 
ing to different ionization channels and their dependence 
on laser intensity. In Figure we have plotted the branch- 
ing ratios of partial ionization yields, namely the number 
of ions left in a particular state of the Magnesium ion 
(ground or excited) . The branching ratios have been calcu- 
lated by integrating over the energy the partial photoelec- 
tron spectrum associated with that state. It is well-known 
that, within lowest-order perturbation theory, the branch- 
ing ratios are independent of laser intensity, as we can ob- 
serve at intensities less than 2 x 10 12 W/cm 2 for 3p and 4s 
ionization thresholds. At moderate intensities, the number 
of ions left in Mg + (3p) state, represented by the dotted 
curve, is greater than the ions left in the Mg + (3d) state 
(dot-dashed curve) and Mg + (4s) state (dashed curve). For 
higher intensities, in the saturation domain, the amount 
of ions left in the Mg + (3p) and Mg + (3<i) excited states 
are of comparable value. 

We also present in Figure [HI the partial photoelectron 
energy spectra corresponding to ions left in the ground 
state (full curve) and in different excited states: 3p (dot- 
ted curve), 3d (dashed curve), 4s (dot-dashed curve) cal- 
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culated for a peak intensity of the laser pulse 2 x 10 12 
W/cm 2 . The spectrum associated with the 3p threshold is 
hidden under the dominant contribution associated with 
the 3s threshold. This overlap can be attributed to the fact 
that the 3s and 3p thresholds are about 4.44 eV apart, so 
that the first photoelectron peak above the 3s threshold 
has only 0.18 eV energy less than the first peak above the 
3p threshold. 

5 Conclusion 

In conclusion, we have calculated ionization yields and 
ATI spectra of Magnesium subject to a pulse of dura- 
tion of 93 fs and 120 fs for the Ti:sapphire laser (804 nm) 
for the range of intensities between 10 12 W/cm 2 and 5x 
10 13 W/cm 2 . For this, we have employed a well estab- 
lished method, based on a configuration interaction dis- 
cretized two-electron states, developed over the last few 
years [H[8 , 20 , 23 , 24] , and in addition we have presented an 
alternative method to solve nonpertubartively the time- 
dependent Schrodinger equation that handles ionization 
stages of two-electrons atoms with multiple continua (many 
open channels). Partial ionization yields, photoelectron 
spectra and photoelectron angular distributions (though 
not presented in this work) may be obtained with this 
open-channel method. 

As already mentioned in the text, the slopes of the 
ion yield as a function of laser intensity that we obtain 
are considerably lower than those observed by Gillen et 
al [5]. Our calculations have not included averaging over 
the interaction volume. Such averaging in any case tends 
to lower rather than increase the slopes of the ion yield 
curves. The effect of resonance with intermediate atomic 
states is inherent in our calculations. But in any case, 
it is not expected to raise slopes. Thus at this point, we 
have no clue as to the possible reasons for the discrepancy. 
Hopefully, further experimental and possibly theoretical 
work will shed some light on this puzzle. 

The work of G.B-Z. was supported by the European Research 
network Program contract No HPRN-CT-1999-00129. One of 
the authors (G.B-Z.) is indebted to Dr. Takashi Nakajima for 
useful discussions. 
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TABLE - FIGURE CAPTIONS 

FIG. 1. Energy-level diagram for Magnesium showir 
the levels relevant for this study. The energies are re 
erenced with respect to the neutral Magnesium grour 
state. 

FIG. 2. Ionization yield in Magnesium as a functic 
of peak laser intensity for photon energies uj = 1.54 e^ 
Pulse duration is 93 fs (35 optical cycles) and 120 fs (4 
optical cycles). 

FIG. 3. Population of 3s 2 , 3s3p, 3s4/ states and ioi 
ization yield (Y), at the end of the pulse, as a function i 
the peak laser intensity, for the 93 fs pulse. 

FIG. 4. Population of 3s 2 , 3s3p, 3s4/ states and ioi 
ization yield (Y) as a function of the interaction time, f< 
the 93 fs pulse, and peak intensity 1.3 x 10 13 W/cm 2 . 

FIG. 5. Comparison of ionization yields in Magnesiu: 
calculated with FRBC method (full curve) and FXB 
(dotted curve) method as a function of the peak laser it 
tensity for photon energy at 804 nm (1.54 eV). 

FIG. 6. Ionization yield in Magnesium as a functic 
of peak laser intensity for photon energies to = 1.54 c^ 
Pulse duration is 93 fs (35 optical cycles). 

FIG. 7. Branching ratios of Mg + ions left in differei 
excited states as a function of laser intensity. The paran 
eters are the same as in Figure O 

FIG. 8. Photoelectron ATI energy spectra correspom 
ing to Mg + left in the ground state as well as a few excite 
states. The laser peak intensity is in = 2 x 10 12 W/cm 
the rest of the parameters are the same as in Figure El 
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